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hole or neutron star mergers, from core-collapse supernovae and even from the in¬ 
flation period in cosmology. These waves exhibit a unique effect called memory or 
Christodoulou effect, which in a detector like LIGO or LISA shows as a permanent 
displacement of test masses and in radio telescopes like NANOGrav as a change in 
the frequency of pulsars’ pulses. It was shown that electromagnetic helds and neu¬ 
trino radiation enlarge the memory. Recently it has been understood that the two 
types of memory addressed in the literature as ‘linear’ and ‘nonlinear’ are in fact two 
different phenomena. The former is due to helds that do not and the latter is due to 
helds that do reach null inhnity. 
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1 Introduction 

Gravitational waves occur in spacetimes with interesting structures from the points 
of view of mathematics and physics. These waves are huctuations of the spacetime 
curvature. Recent experiments [Ij claim to have found imprints in the cosmic mi¬ 
crowave background of gravitational waves produced during the inhation period in the 
early Universe, though it has been argued |2H] that foreground emission could have 
produced the experimental signal. Gravitational radiation from non-cosmological 
sources such as mergers of binary black holes or binary neutron stars as well as core¬ 
collapse supernovae are yet to be directly detected. We live on the verge of detection 
of these waves and thus a new era of astrophysics where information from formerly 
opaque regions will be laid open to us. Experiments like LIGO on Earth or LISA 
planned for space or the pulsar-timing arrays using radioastronomy like NANOGrav 
are expected to detect gravitational waves in the near future. These waves exhibit 
a memory effect, called the Ghristodoulou effect, which in the LIGO or LISA ex¬ 
periments manifests itself by displacing test masses permanently and in the radio 
telescopes like NANOGrav by changing the frequency of pulsars’ pulses. The follow¬ 
ing works by Zel’dovich, Polnarev, Braginsky, Grishchuk and Thorne m, n. ra 
treat a ‘linear’ memory due to a change in the second time derivative of the source’s 
quadrupole moment. This effect was believed to be small. However, Ghristodoulou 
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ini established the ‘nonlinear’ memory effect resulting from energy radiated in gravi¬ 
tational radiation. Christodoulou also showed that his new effect is much larger than 
the formerly known ‘linear’ one. The hrst two of the present authors showed 0. i 
that, what in the literature has been known as the ‘linear’ memory is due to helds 
that do not reach null inhnity, and what has been referred to as ‘nonlinear’ memory 
is due to helds that do reach null inhnity. Christodoulou irzi established his result 
for the Einstein vacuum equations. In general relativity (GR) matter or energy helds 
on the right hand side of the Einstein equations Q generate curvature. Do these 
helds contribute to the memory ehect? And how large would such a contribution 
be? The hrst and the third of the present authors with Chen proved 0 . 0 that 
electromagnetic helds enlarge the Christodoulou ehect. The hrst and the second of 
the present authors established [8] a contribution to memory from neutrino radiation 
in GR, where the neutrinos are described by a null huid. The same authors also 
investigated memory for more general energy-momentum tensors |H], and they found 
an analog to gravitational memory within the regime of the Maxwell equations, 
which establishes this phenomenon for the hrst time outside of GR. 

Radiative spacetimes are solutions of the Einstein equations 

Rfiu - ( 1 ) 

with denoting the Ricci curvature tensor, R the scalar curvature, the metric 
tensor and the energy-momentum tensor of any helds considered. If the right 
hand side is not trivial, then corresponding matter equations have to be provided. 
Denote by (M, g) our 4-dimensionaI spacetime manifold where g is a Lorentzian 
metric obeying Q. 

Often the left hand side of Q is denoted by which is called the Einstein 
tensor. The twice contracted Bianchi identities yield 

D‘'G,y = 0 . ( 2 ) 

We see from equations Q and (|^ that the energy-momentum tensor is symmetric 
and divergence free 

= 0. (3) 

The Einstein equations ([^ are of hyperbolic character and exhibit interesting 

geometric structures. In mathematical GR the methods of geometric analysis have 
proven to be most ehective in solving some of the most challenging burning problems. 
It turns out that a lot of information about radiation can be gained from the study 
of asymptotics of spacetimes that were constructed within stability proofs. Whereas 
in the latter, the initial data has to fulhll corresponding smallness assumptions, the 
emerging asymptotic structures are independent of smallness. In fact, these results 
hold for large data such as binary black hole mergers and other sources of gravitational 
radiation. We address this connection in sections |3] and HI 

Many physical cases require to study the Einstein equations in vacuum. By setting 
to zero the right hand side of Q we obtain the Einstein vacuum (EV) equations 

Rm- = 0 . (4) 
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Solutions of the EV equations are spacetimes {M,g), where M is a four-dimensional, 
oriented, time-oriented, differentiable manifold and g is a Lorentzian metric obeying 
the EV equations. This is the hrst of the following matter models. Among the most 
popular ones we hnd: 

1. The Einstein vacuum (EV) case with = 0 and equations (|^. 

2. The Einstein-Maxwell (EM) case with where 

F denotes the electromagnetic held, thus an antisymmetric covariant 2-tensor. 
Tpy being trace-free, the Einstein-Maxwell system reads 

Rpu = 

D^Fp. = 0 (5) 

D^^*Fp, = 0 . 

3. The Einstein-Yang-Mills (EYM) equations with Lie group G with Tp^ being the 
stress-energy tensor associated to the Yang-Mills curvature 2-form Fp^. The 
latter takes values in the Lie algebra of G and has to satisfy the EYM system. 
F is written in terms of the potential 1-form A. The tensor Tp^ is trace-free as 
in the Maxwell case. But now for non-Abelian groups the YM potential does 
not disappear from the equations. 

4. Fluid models such as a perfect huid, in particular a null huid. In each case 
Tpy takes a specihc form. The obvious conservation laws have to be fulhlled 
including the particle conservation law. 

Section gives a brief history of the efforts to detect gravitational waves. In 
section]^ we introduce the Cauchy problem for the Einstein equations and explain the 
stability results for asymptotically flat spacetimes. Section connects these results 
to gravitational radiation. In part we investigate the Christodoulou memory effect 
of gravitational waves for the Einstein equations. Finally, in the last two sections 
and the prospects for the detection of gravitational wave memory is discussed and 
a brief discussion of gravitational waves in the cosmological setting is given. 

As the Cauchy problem, black holes, stability of the Kerr solution, cosmology and 
related topics are covered by other articles, we only briefly touch these here and say 
as much as needed for our purposes. Also, for a complete list of references regarding 
these topics we refer to the corresponding articles in this volume. 

2 History of Gravitational Wave Observations 

Gravitational radiation can be understood in analogy with the more well known 
case of electromagnetic radiation. The motion of electric charges makes waves of 
changing electric and magnetic helds which propagate outward from their sources, 
travel at the speed of light, carry energy away from their sources, and fall off like 
r~^ at large distances r from the source. Electromagnetic waves can be detected 
directly through their effect on the motion of charges, and indirectly by measuring 
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the amount of energy that their sources lose. Similarly, the motion of masses makes 
waves of changing spacetime curvature which propagate outward from their sources, 
travel at the speed of light, carry energy away from their sources, and fall off like r~^ 
at large distances r from the source. Gravitational waves can be detected directly 
through their effect on the motion of objects in free fall, and indirectly by measuring 
the amount of energy that their sources lose. The main differences between the two 
types of radiation are (1) gravitational waves are much weaker, and (2) the equation 
that describes gravitational radiation is nonlinear. Because of the weakness of the 
gravitational interaction, there are no detectable manmade sources of gravitational 
radiation: any gravity waves that we could produce in the lab would be too small 
to be detected. Only astrophysical sources have any hope of producing detectable 
gravitational waves, and the most promising sources are those in which gravity is 
strongest and most dynamical: colliding black holes and neutron stars, and supernova 
explosions. 

Despite the crucial nature of nonlinearity in this chapter, we will hrst consider 
what happens when we neglect it. That is we will consider that class of solutions 
where the gravitational helds are everywhere sufficiently weak that they can be well 
approximated by neglecting the nonlinear terms in the Einstein held equation. Then 
the held equation reduces to a linear equation that is similar to Maxwell’s equations 
for electromagnetic helds. In addition we further restrict attention to the case where 
the sources are moving slowly compared to the speed of light. Then as for Maxwell’s 
equations one obtains an expansion for the radiation held in terms of time derivatives 
of the multipole moments of the source. [2^ For electromagnetism the conservation 
of charge forces the hrst term in this expansion to vanish and so the dominant term 
involves the dipole moment. For gravity the conservation of energy and momentum 
forces the hrst two terms in the expansion to vanish, and the dominant term involves 
the quadrupole moment. In particular, in this weak held slow motion approximation 
for gravity one obtains the following formula for the rate at which energy is lost in 
gravitational radiation: 

1 * k* ah * k,* 

P=lQ Qal ( 6 ) 

Here P is the power radiated in gravitational radiation, = J d^xpx°‘x^ is the mass 
quadrupole moment tensor, is the trace-free part of and an overdot denotes 
derivative with respect to time. This is the so called quadrupole formula of general 
relativity. As derived within linearized gravity, the quadrupole formula applies only 
to systems whose motions are caused by non-gravitational forces. However, there is 
also a derivation of the quadrupole formula using the post-Newtonian approximation 
(see e.g. 0) that applies also to systems whose motions are inhuenced by gravity. 
The quadrupole formula immediately suggests a strategy for the indirect detection 
of gravitational waves: hnd an astrophysical system whose motion can be measured 
very accurately. Measurement of the motion allows one to calculate the quadrupole 
moment and thus the rate at which general relativity predicts that the system will 
lose energy. But that loss of energy will give rise to a change in the motion of the 
system, which can be measured through prolonged observation of the system. It is 
precisely this strategy that was adopted by Hulse and Taylor [US] in their discovery and 
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observations of the binary pulsar. This system consists of two neutron stars in orbit 
around each other. One of these neutron stars is a pulsar: a rotating neutron star 
from which we receive a pulse of radio waves once each rotation. This series of pulses 
acts as an extremely accurate clock which also enables measurement of all details 
of the orbit through Doppler shifts in pulse arrival times from the orbital motion 
of the neutron star. As the binary pulsar system loses energy, the size of the orbit 
shrinks and the period of the orbit gets shorter. The quadrupole formula for energy 
loss then yields a prediction for the rate of change of the period. Hulse and Taylor 
measured the rate of change of period and found agreement with the prediction of 
general relativity. In so doing they have (indirectly) detected gravitational radiation. 

Though the work of [66] leaves no doubt about the existence of gravitational 
radiation, one would still like to directly detect it. To understand the strategies for 
direct detection, it is helpful to turn to the properties of geodesics. Objects in free 
fall follow spacetime geodesics. The spatial separation between two nearby geodesics 
will change with time due to the nonuniformity in the gravitational held encoded in 
the spacetime curvature. Specihcally the separation of geodesics is governed by the 
geodesic deviation equation 

y = -RtujS^ (7) 

Here one geodesic dehnes the rest frame that determines the time coordinate t and 
spatial coordinates x\ while S'* is the spatial separation between the two geodesics 
and Rafj'yS is the Riemann tensor. It is instructive to look at the behavior of eqn. ([^ 
in the weak held slow motion case. We hnd 

[Qy] (8) 


Here V\\ denotes ‘projected orthogonal to the radial direction and trace-free.’ While 
a gravitational wave changes the positions of objects in free-fall, one might have 
expected that after the wave has passed the objects return to their original positions. 
However, using eqn. ([^ in eqn. ([^ and integrating twice with respect to time one 
obtains 


AS’ 


= -S^V 

r 


AQ 




(9) 


Here A denotes the diherence in a quantity between the begining and end of the 
time of the passage of the gravitational wave. Thus there is a change in the positions 
of the free fall objects: a gravitational wave memory. Note however, that eqn. 
only holds for slowly moving sources and in the linearized approximation to general 
relativity. For this reason, Christodoulou calls this contribution to gravitational wave 
memory the “linear memory.” When one does not make the linear approximation 
and instead uses the full nonlinear theory of general relativity there is an additional 
contribution to the memory. 

The geodesic deviation equation immediately suggests two strategies for direct de¬ 
tection of gravitational radiation which we will call “bar” and “beam.” The geodesic 
deviation equation implies that the effect of a gravitational wave on a somewhat rigid 
object (a “bar”) will be to make the bar vibrate. This effect will be especially strong 
if the gravitational wave has a frequency that is close to the resonant frequency of the 
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bar. One then isolates the bar from all extraneous sources of vibration (to the extent 
possible) and measures the tiny vibration induced by the passage of gravitational 
waves. This is the strategy used by Weber [72] in his attempt to detect gravitational 
waves. Unfortunately, this strategy was not successful, and in hindsight it is easy 
to see a major weakness of the bar strategy: since the right hand side of eqn. ([^ 
contains , it follows that the signal (change in separation) is proportional to the 
original separation. Thus, to get a signal that is sufficiently large to measure, one 
should start with geodesics that have a sufficiently large separation. There is a prac¬ 
tical limit to the size of the bar and thus a practical limit to the sensitivity of the 
bar detector. 

This difficulty immediately suggests the alternative strategy: have the detector 
consist of widely separated objects and keep track of their spatial separation using 
light (“beams”). One especially accurate way of using light to measure spatial sepa¬ 
ration is laser interferometry: a laser interferometer consists of a laser, two mirrors, 
and a half reflective mirror called a beam splitter, with the mirrors and beam splitter 
placed in an L shaped conhguration. The laser shines on the beam splitter which 
splits the beam into two pieces, each of which travels down one arm of the L, reflects 
off the mirror at the end of the arm and goes back to the beam splitter where the 
two halves of the beam recombine. Since the two arms of the L are not exactly equal 
length, the distances traveled by the two light beams are different. If this path differ¬ 
ence is an odd number of half wavelengths of the laser light then the two light beams 
will interfere with each other (destructive interference) and the combined beam will 
be completely dark. A laser interferometer set at destructive interference allows a 
very sensitive measurement of any change in the position of the mirrors, since any 
change in position of a fraction of a wavelength of light will lead to the dark being 
replaced by a partial beam, and from the brightness of the partial beam one can read 
off the fraction of the wavelength that the mirrors have moved. This beam strategy 
is used by LIGO (Laser Interferometer Gravitational wave Observatory). [53] LIGO 
consists of two detectors (One in Hanford, WA and one in Livingston, LA), each a 
laser interferometer with a 4 km arm length. LIGO began operations in 2001 and has 
not so far detected any gravitational waves. However, the detector has recently been 
upgraded (more powerful laser, better mirrors, better shielding of the mirrors from 
extraneous vibration, etc.). With the upgrades one can fairly conhdently expect the 
hrst direct detection of gravitational waves within the next few years. Other laser 
interferometers include VIRGO (located in Gascina, Italy with 3 km arm length) and 
GEO600 (located in Sarstedt, Germany with 600 m arm length). 

If bigger interferometer arm length is better, then one can get it by going to space. 
Here the strategy is to use laser interferometry, but with the lasers and mirrors 
mounted on three satelites in a triangular conhguration. This is the idea behind 
the proposed LISA (Laser Interferometer Space Antenna) project. LISA has been 
variously proposed as a NASA project, an ESA (European Space Agency) project, 
and as a joint project between NASA and ESA. It may eventually be built and 
launched, but not any time in the near future. Even larger “arm length” can be 
obtained through pulsar timing. When a gravitational wave passes between a pulsar 
and the Earth, it distorts spacetime and affects the timing of the arrival of the pulses 
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at the Earth. [23] If several pulsars are being timed, then the gravity wave will affect 
the timing of all their pulses, and by an amount that depends on the direction of 
the line of sight to the pulsar as well as on the direction and the polarization of the 
gravitational wave. Thus, the signature of a gravitational wave is a change in the 
timing of several pulsars with a particular angular pattern. This is the idea behind 
the NANOGrav (North American Nanohertz Observatory for Gravitational waves) 
project. Radio astronomers have done pulsar timing for many years, but recently 
there has been a concerted effort to use pulsar timing to detect gravitational waves. 
Here the more pulsars that are timed and the longer they are observed, the greater 
is the overall sensitivity to a gravitational wave signal. 

3 The Cauchy Problem and Stability 

The Gauchy problem for the Einstein equations is discussed first. Then stability 
results are explored. 

3.1 The Cauchy Problem 

The Einstein equations split into a set of constraint equations which the initial data 
has to satisfy and a set of evolution equations. Depending on the matter and energy 
present on the right hand side of ([^ we also specify corresponding equations for 
these helds. The goal is to construct a spacetime by solving an evolution problem 
of the Einstein equations. An initial data set consists of a 3-dimensional manifold 
H, a complete Riemannian metric g, a symmetric 2-tensor k and a well specified 
set of initial conditions corresponding to the matter-fields. These have to satisfy the 
constraint equations. 

A Cauchy development of an initial data set is a globally hyperbolic spacetime 
(M, g) satisfying the Einstein equations together with an imbedding i : H ^ M such 
that A(^) and A(A;) are the first and second fundamental forms of i{H) in M. 

Most commonly used are two settings of foliations. In the hrst setting, the space- 
time is foliated by a maximal time function t and an optical function u, respectively. 
The time function t foliates our 4-dimensional spacetime into 3-dimensional space- 
like hypersurfaces Ht, being complete Riemannian manifolds. Whereas the optical 
function u induces a foliation of {M,g) into null hypersurfaces C^, which we shall 
refer to as null cones. The intersections Ht f] Cu = St^u are 2-dimensional compact 
Riemannian manifolds. In the second setting, the spacetime is foliated by an optical 
function u and by a conjugate optical function u. The resulting outgoing and incom¬ 
ing null hypersurfaces constitute the double-null foliation. The characteristic initial 
value problem then consists of initial data given on an initial null cone instead of a 
spacelike hypersurface and the corresponding evolution equations. In this paper, we 
write the equations with respect to the hrst setting. 


7 


The evolution equations in the EV situation read: 


^9ij 

dkij 


2^kij 


ViVj<h - {Rij + kij trk - 2kimk’f)^ 


whereas the constraint equations are given as: 


— Vj trk = 0 

R + {trkY — \k\^ = 0 . 

Here, <h denotes the lapse function $ := (—g'^Witdjt)~^, and V is the covariant 
derivative on H. Moreover, barred curvature quantities denote the corresponding 
curvature tensors in H. If the time function t is maximal, then trk = 0. The time 
vector held T, thus the future-directed normal to the foliation is given by T* = 
— t^'^g'-tdjt and Tt = T^dit = 1. The components of the inverse metric are g'’^ = 
{g~^y^. The metric reads g = + g. The foliations by t and u supply us 

with natural vectorhelds to work with. Consider St^u that we introduce above as the 
intersection of Ht and Cu- Let N be the spacelike unit normal vector of St^u in Ht- 
Let {ea}a=i ,2 be an orthonormal frame on St^u- Then {T, N,e 2 ,ei) is an orthogonal 
frame. Also, let T = denote the unit future-directed time vector, then consider 

the pair of null normal vectors to St^u, namely L = T + N and L = T — N. Together 
with {ea}a=i, 2 , they form a null frame. 

In this article, we consider asymptotically hat spacetimes in GR. These describe 
isolated gravitating systems such as galaxies, neutron stars or black holes. These 
are solutions where (M, g) approaches hat Minkowski space with diagonal metric 
rj = (—in the limit of large distance from a spatially compact region. 
Thus, they can be thought of as having an asymptotically hat region outside the 
support of the matter. 

Consider an asymptotically flat initial data set, i.e. outside a sufficiently large 
compact set K, H\K is diheomorphic to the complement of a closed ball in and 
admits a system of coordinates in which g —Sij and A: —?• 0 fast enough. 

On the way to the well-posedness of the Cauchy problem for the Einstein equa¬ 
tions, we hrst cite de Bonder’s [22] work on harmonic gauge. Several authors then 
contributed in various ways such as Friedrichs, Schauder, Sobolev, Petrovsky, Leray 
and others, which led to the existence and uniqueness theorems for general quasilin- 
ear wave equations. In [39] Leray formulated global hyperbolicity. Finally, the hrst 
proper formulation of the Cauchy problem for the Einstein vacuum equations is due 
to Y. Choquet-Bruhat. In 1952, Y. Choquet-Bruhat [13] treated the Cauchy prob¬ 
lem for the Einstein equations locally in time establishing existence and uniqueness 
of solutions, reducing the Einstein equations to wave equations, introducing wave 
coordinates. The local result led to a global theorem proved by Y. Choquet-Bruhat 
and R. Geroch in [15], stating the existence of a unique maximal future development 
for each given initial data set. As a logical further question arises: Is this maxi¬ 
mal future development complete? No, not necessarily. This answer by R. Penrose 




was formulated in his incompleteness theorem. The latter says that ‘generic’ initial 
data containing a closed trapped surface will yield a spacetime M which is not null 
geodesically complete (towards the future), in particular it will have a maximal fu¬ 
ture development which is incomplete. A closed trapped surface S' in a non-compact 
Cauchy hypersurface if is a two-dimensional surface in if, bounding a compact do¬ 
main such that 

trx <0 on S . 

X is the second fundamental form of a spacelike surface S with respect to the outgoing 
null hypersurface C properly dehned in section The theorem of Penrose and its 
extensions by S. Hawking and R. Penrose led to ask the following: Is there any non¬ 
trivial asymptotically flat initial data whose maximal development is complete? Yes. 
This answer was given by D. Christodoulou and S. Klainerman [20] and is discussed 
in the next section. Instead here, we briefly turn to the opposite direction of black 
hole formation. 

In 2008 D. Christodoulou [I9| proved that by the focussing of gravitational waves 
a closed trapped surface and therefore a black hole will form for initial data that 
does not contain any closed trapped surface. Since then several authors including S. 
Klainerman, J. Luk, I. Rodnianski, P. Yu, contributed to this line of research. See 
the corresponding articles on this topic for detailed information and references, the 
most recent being [33] . 

Recent regularity results by S. Klainerman, I. Rodnianski and J. Szeftel [|36j. 133. 
[38] . [62] . [63] , [M] , [65] ] proved the bounded curvature conjecture. In particular, 
the authors showed for solutions of the EV equations that the time of existence 
depends only on the L^-norm of the curvature and a lower bound of the volume 
radius of the corresponding initial data set. 

3.2 Stability 

In proving the global nonlinear stability of Minkowski space [20] D. Christodoulou and 
S. Klainerman achieved the global result. Later, a proof under stronger conditions 
for the global stability of Minkowski space for the EV equations and asymptotically 
flat Schwarzschild initial data was given by H. Lindblad and I. Rodnianski 03. 03. 
the latter for EV (scalar held) equations. They worked with a wave coordinate 
gauge, showing the wave coordinates to be stable globally. Concerning the asymptotic 
behavior, the results are less precise than the ones of Christodoulou-Klainerman [CK] 
in m- Moreover, there are more conditions to be imposed on the data than in [20] . 
There is a variant for the exterior part of the proof from [20] using a double-null 
foliation by S. Klainerman and F. Nicolo in [Mj. Also a semiglobal result was given 
by H. Friedrich [22] with initial data on a spacelike hyperboloid. N. Zipser eg. eg 
generalizes the CK proof [2D] to the Einstein-Maxwell case. L. Bieri 0.131 generalizes 
the CK results m for the Einstein vacuum equations by assuming only 3 derivatives 
of the metric to be controlled as opposed to 4 in [2Q] and by assuming that the initial 
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data as r —>■ cxD behaves like 


gij = Sij + 03{r 2 ) 

kij = 03(r“5) (10) 

which relaxes the decay of 1201 by one power of r. The latter establishes the border¬ 
line decay of the data. We recall that Christodonlon-Klainerman in [20] considered 
strongly asymptotically flat initial data which as r —)■ cx) take the form 

2M _3 

Qij = (1 + —) dij + 04 (r 2 ) 

kij = (11) 

where M denotes the mass. 

Farther resnlts in varions directions were obtained in la, 031, 01 , 051, and 
in the cosmological context [55], [6T], [60], [56], [16]. For more details the reader 
may consult the papers on the corresponding topics in the present volume. In this 
article we concentrate on the stability results and methods which are important for 
radiation. 


Theorem 1 [CK in Strongly asymptotically flat initial data {HQ,gQ,kQ) of the 
type 0 which is sufficiently small, yields a unigue, causally geodesically complete, 
globally hyperbolic solution {M,g) of the EV eguations. The spacetime {M,g) is 
globally asymptotically flat. 


The proof of this celebrated result as well as the proofs of [75], [TB] and 0.0 use 
geometric-analytic methods, thus are invariant, in particular they do not depend on 
any coordinates. All these proofs use the spacetime-foliation introduced above in 
section 3T as the hrst setting, where the foliations are induced by a maximal time 
function t and an optical function u. These proofs rely on ‘the right energies’ which 
we hnd in the Bel-Robinson tensor. Very roughly the proofs can be summarized in 
the following steps: 1) Introduce and estimate the energies, that can be controlled. 2) 
Estimate the curvature by a comparison argument with these energies. 3) Estimate 
all the geometric quantities within a bootstrap argument by curvature assumptions 
and evolution equations. Make use of Hodge systems on St^u coupled to transport 
equations along the null and spacelike hyper surf aces. 

Note that in the situations of BD], ivg, 051 one sees radiation, whereas in the 
situation of 0.0 radiation cannot be seen due to the lower decay of the data. 

The null asymptotic picture of these results is independent from the smallness 
assumption in the theorems. Whereas smallness guarantees existence and unique¬ 
ness of solutions, the asymptotic results hold true for large data as well. Thus, the 
asymptotics for typical sources of gravitational waves such as merging black holes 
generating large data can be studied using these results. 

In the next section, we investigate the null asymptotics of radiative spacetimes. 
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4 Radiative Spacetimes and Null Infinity 

In order for a spacetime to carry radiation, it cannot be spherically symmetric. In 
what follows, we consider spacetimes with asymptotic structures as worked out by 
Christodoulou and Klainerman in IZDl, by Christodoulou in DU for the EV case, as 
by Zipser in [75], [76], by the hrst and third present authors with Chen in 0.0 in 
the EM situation and as by the hrst and second present authors in [S], and by the 
hrst author in [ 3 ] in the neutrino case. 

The most important quantities to describe radiation are the shears of incoming 
respectively outgoing null hypersurfaces. These shears x and x are the traceless parts 
of the corresponding second fundamental forms. The second fundamental form x of a 
spacelike surface S relative to the outgoing null hypersurface C is given by x(X, V) = 
g{DxL, Y) for any X, Y in TpS and L generating vectorheld of C. Correspondingly, y 
is dehned in the analogous way with L being replaced by L the generating vectorheld 
of the incoming null hypersurface C. 

These shears have limits along each null hypersurface as t —?• oo. In all the 
works cited in the hrst paragraph of this section it is proven for the corresponding 
situations that 

lim r^x = S(M, •), ^hm ry = H(m, ■) . ( 12 ) 

Versions of the following theorem in radiative spacetimes are proven in the works 

izni. na, ca, ng, 0 . 0 , 0 , a. 

Theorem 2 The limit H(m, •) behaves like 

l‘='( m,^ (1 + |m|) 

where \ ■ |o denotes the pointwise norm on with respect to the metric 7 being the 
limit as t ^ 00 for each u of the induced metrics on St^u rescaled by ^. 

The behavior of E(m, •) is given by Theorems 2, 3 and 4. In particular, S tends to 
limits S’*" and with u —)■ +cxd and u —)■ —00 respectively, which is guaranteed 
by Theorems 2 and 3. Moreover, it follows from Theorem 4 that these limits and in 
particular the diherence — S“) are nonzero. 

With respect to the (t, u) foliation introduced above we dehne the following cur¬ 
vature components, with index notation A,B^ { 1 ) 2 }, 

aAB = R{eA, T, cb, L) , = R^ca, T, cb^L) 

13a =-R{eA, L, L, L) , =-R{eA, L, L, L) (13) 

p=^R{L,L,L,L) , a=^*R{L,L,L,L) 

In the radiative spacetimes of interest all these components except for oab and /3a 
converge to their well-dehned limits along as t —)■ cxd for each u. The exceptions 
aAB and /3a tend to zero. The most important limit in view of the following discussion 

is 

Awiu, •) = lim ra (14) 

Cu —^00 
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where the index W denotes the Weyl curvature component to be distinguished from 
curvature generated through matter helds in the Einstein equations. When working 
with the EV equations (|^ the Weyl curvature is equal to the Riemann curvature. 
However, as soon as helds are present on the right hand side of the full Einstein equa¬ 
tions Q they manifest themselves in the energy-momentum tensor generating 
extra curvature. Recall that the Weyl tensor (which is traceless) is computed from 
the Riemann-, Ricci- and scalar curvature tensors as 

1Rq:/37(5 Rafi'^S “h QfisRa'y 9/3'yRaS dasRlB'y') R g (^Qa-yQjSd 9a59/3'y') R ■ 

Denote by the component of the energy-momentum tensor and its null limit 
by 7cn = limc^ For the EM equations Tc^ is a quadratic of the electro¬ 

magnetic held component which decays like r~^ plus extra decay in u. In the case 
of a null huid Tc^ again exhibits the decay behavior of plus decay in u. In all 
the corresponding cases investigated in ra, mi, iTg, izsi, 0, i, 0, a there is an 
interesting relation between the limiting shears and the Weyl curvature limit: 


Theorem 3 The curvature component Aw o,nd the shears E, S are related through 

BE __l 
du 4 

Moreover, it is \Awiu-, •)! 

By 4' denote a Ricci rotation coefficient of the null frame, including trx, trx, the 
shears y, x ^^nd the torsion 1-form (a = \g{,DiJj, ca) plus its conjugate. By $ denote 
a component of the Weyl curvature and by T a component of the energy-momentum 
tensor. fjAB is the traceless tensor component of the second fundamental form k 
tangential to the Rt^^-surface. ^ is the covariant derivative on St^u- Then it is 


f ^X = + “4'” + “T” 

f = “TT” + “T(8)T” + + “T” 


Propagation equations for trx and fry: 


dtrx 

ds 

dtrx 

ds 

The Gauss equation reads 


-itrytry-2p + 2|C|^ 

{Rxf - \xf -87 tT. 


K 


-^trxtrx + ^X-X - ‘W” + 


contributions from T 


(16) 

(16) 
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The null Codazzi and conjugate null Codazzi equations read 


dijvx = -X-C++ - ‘W” 

+ contributions from T 

dijvx = X - C+]^{ytrx - (trx) + ‘W” 

+ contributions from T 

These equations together with other structure equations, which we do not cite, are 
used to derive radiation laws. To start the investigation, we dehne the Hawking mass 

r r f 

m{t,u) = 2 + ^ / ^^Xtrx- 

By taking ^ and ^ of m and investigating the null limits, we derive the Bondi mass 
and the Bondi mass loss relation at null inhnity. First, we hnd that the the Hawking 
mass m{t,u) tends to a limit M{u) called the Bondi mass along each Cu as t ^ oo. 
Second, Bondi mass loss formulas are derived. In the EV case ra, HZ] this takes the 
form 



Whereas in all the other cases [zg, [zi, IS], 0,0, in this formula exhibits an extra 
term 

The positive constant C and the structure of the positive limit 7c„ depend on the 
specihc situation and are investigated in the works cited. 

Dehning the function 

/ OO 

\E\^ + CrcJu (17) 

■OO 

T/Stt is the total energy radiated away per unit solid angle in a given direction. This 
seeds the basis for section IH 

4.1 Peeling Behavior 

To understand radiating properties of a spacetime, we need to know the behavior of 
the curvature at null inhnity. Thus decomposing this tensor according to a null frame 
allows one to describe this behavior. However, what is the asymptotic structure? 

Nowadays, one likes to group the answer to this question into the Christodoulou- 
Klainerman picture and the Newman-Penrose picture. 

The null asymptotic behavior of the Riemann curvature tensor RafSyS is crucial in 
the present discussion. Its specihc decay behavior known as “peeling estimates” was 
discussed by Bondi et ah m, Sachs [57| and Penrose (SO] . In their pioneering paper 
|1H] Newman and Penrose use a null frame to decompose the important geometric 
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quantities. Even earlier [5T] Pirani finds peeling properties for curvature compo¬ 
nents. These are important for spacetimes which can be conformally compactihed 
(see subsection 4.1.1). 

In [20] Christodoulou-Klainerman derive very precise estimates for the Riemann 
curvature, they obtain peeling up to one component. See Christodoulou’s paper 
[16] for further details on these spacetimes. (See subsection 4.1.2] ) However, global 
solutions obtained by Klainerman-Nicolo [3l] when assuming stronger decay and 
regularity were shown to possess peeling estimates [33] . 

The Newman-Penrose formalism refers to the structure equations of a null frame 
coupled with the Bianchi identities. 


4.1.1 Peeling in the Newman-Penrose Picture 


Peeling is addressed in the following pioneering papers [IH], [H], [50] . [51], [57]. If 
one imposes a smooth conformal compactihcation of the spacetime, then a specihc 
hierarchy of decay follows for the curvature components with respect to the frame 
considered. In fact, “smooth” means at least We can think of this conformal 
compactihcation as introducing a boundary future null inhnity. Based on this one 
can dehne rescaled Newman-Penrose curvature scalars being functions on this bound¬ 
ary. These are related to rescaled spin coefficients. In fact, often in literature the 
terminology “spin-coefficient formalism” is used for “Newman-Penrose formalism”. 


Given such a conformal compactihcation, the curvature components (13) show 
the following decay in r as f —)■ cx) for hxed u\ 


a = 0{r , P = 0{r ^), 

p,a = 0(r“^) , 

[3 = , a = 


4.1.2 Peeling in the Christodoulou-Klainerman Picture 

In [20] Christodoulou-Klainerman derive the following behavior in r for the curvature 
components (13) as t —)■ cxd for hxed u: 

a = 0(r~^) , /3 = 

p,a = 0(r-3) , 

(3 = 0(r“5) ^ a = 0(r“5). 


We note that in [SU] the components a and f3 of the curvature fail peeling. In fact, 
the results of [20] do not allow for a compact but for a complete future null inhnity. 
This behavior is a direct consequence of the stability theorem [T] Now one could ask 
if f3 were to exhibit the stronger fall-oh 0{r~^), if only one asked for more decay of 
the initial data in theorem Christodoulou [TU] shows that this is not the case for 
generic initial data. He [H] proves that at best what one can achieve is f3 behaving like 
r“^logr. The reason is that the curvature component (3 is coupled to the radiative 
amplitude, that is to the shear y, and by investigating their limits at null inhnity 
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/3 would fulfill peeling only if a certain expression involving the past limit of this 
shear vanishes. However, it can be computed that for generic inital data this very 
expression cannot vanish. 

From this discussion we conclude that the question of peeling can only be un¬ 
derstood from the point of view of a Cauchy problem for the Einstein equations. 
What assumptions are required on the initial data so that a reasonable dehnition of 
future null inhnity follows for the emerging spacetime? This indeed is achieved in [20] 
establishing the latter asymptotics and thereby dehning complete future null inhnity. 


5 ‘Linear’ and ‘Nonlinear’ Memory 


Gravitational radiation, namely the huctuation of curvature, manifests itself in the 
displacements of geodesics. Whenever a binary merger or a core-collapse supernova 
occurs, a gravitational wave packet is sent out which travels at the speed of light along 
outgoing null hypersurfaces. These waves have two effects on geodesics: First the so- 
called instantaneous displacements during the ‘hy by’ of the wave packet, second 
the permanent displacements after passage of the wave packet called the memory. 
This memory has a ‘linear’ and a ‘nonlinear’ part, the latter being the Christodoulou 
memory. 

In the relevant works cited in the previous section, it is proven that 

|S| —)■ 0 as |m| —)■ cxD 

whereas S takes limits S’*" as u —)■ -|-oo and as u —?• —oo. It is also shown that 

E(m) = S- + ^ [ E{u') du' (18) 

^ J — OO 



S(m) — is related to instantaneous displacements and to permanent dis¬ 

placements of geodesics. 

In [20], [I2| for EV, jS], |S| for EM, [S], |1] for Einstein equations with neutrino 
radiation via a null fluid, it is proven that the permanent displacement related to 
E"*" — is governed by the energy S. Generally, these results can be formulated as 
follows. 


Theorem 4 Denote by £ the mean value of the energy £ in (11) on . Denote by 

O _ 

{]) an operator on S^. Let h be the solution with h = 0 on of the equation 


h = £-£. 

Then S’*' — is determined by the following equation on 8“^: 


di/u (E+ - E-) =f h 


( 20 ) 
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Thus, the nonlinear memory of gravitational waves is governed by the energy function 
(Hill- We see that the shear S contributes, and is the only contribution in the EV case 
ini, but also there are contributions through the energy-momentum component Tc^ • 
The contributions from the EM equations are derived in 0 , 0 , for Einstein-null-fluid 
and neutrino radiation in [8] and for a general energy-momentum tensor in [9]. 

Let us now inquire how these results impact geodesic behavior and the spacetime 
geometry. 

We consider three geodesics in spacetime and name them Qq, Qi, ^ 2 - Let T denote 
the unit future-directed tangent vectorheld of Qq and t the arc length along Qq. More¬ 
over, for each t denote by Ht the spacelike, geodesic hyperplane through ^o(^) which 
is orthogonal to T. Then we consider the orthonormal frame held {T, Ei, E 2 , E^) 
along Qq, where {Ei, E 2 , E^) is an orthonormal frame for Hq at ^o(O)) parallely prop¬ 
agated along Qq. This yields for every t an orthonormal frame {Ei, E 2 , E^) for Ht 
at ^o(^)- Next, to a point p lying in a neighbourhood of Qq we assign the cylindri¬ 
cal normal coordinates {t,x^,x‘^,x^), based on Qq, if p E Ht and p = expX with 
X = ^ Under assumptions valid on the Earth, one replaces the 

geodesic equation for Qi and Q 2 by the Jacobi equation (geodesic deviation from Qq). 

M k 

- RkTlT (21) 

with 

RkTlT = R (Ek^T, Ei,T) . 

To see better the structures of the curvature terms, we switch to the null vectorhelds 
again where L = T — E^, L = T + E^. The leading components of the Weyl curvature 
are 


~ R {Ea, L, Eb, Ll) 

_ '^ab{W) ^ ^ _ 2 ^ 

^ABi^) — + o [r ) . 

r 


( 22 ) 

(23) 


Whenever helds are present on the right hand side of the Einstein equations ([^ and 
they have the ‘right’ decay properties we hnd the leading order component of the 
energy-momentum tensor to be 


Tll 


Tc 


Cn 


+ l.o.t. 


(24) 


In order to determine displacements of the geodesics, we have to integrate twice 
equation (21). Here, it is crucial that we know the relations between the two shears 
and the Weyl curvature component Aw- By virtue of theorem]^ and the geometric- 
analytic properties of the relevant spacetimes we derive the permanent displacement 
Ax as 


A a: = - (—) (E+ - S") . (25) 

r 

The right hand side is understood by theorem and equation (0. The case of the 
EV equations was investigated by D. Christodoulou mi, for the EM equations by 
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L. Bieri, P. Chen, S.-T. Yau [5], j^, for the Einstein-null-fluid describing neutrino 
radiation by L. Bieri, D. Garhnkle [8] and for a general energy-momentum tensor by 
the latter authors in [H]. 

6 Detection of Gravitational Wave Memory 

What are the prospects for detection of gravitational wave memory? This question is 
treated in|2n|. Here it is important to remember that there has not yet been any direct 
detection of gravitational waves. Therefore it is best to start by considering whether 
memory is harder to detect than other aspects of gravitational waves. Because of the 
weakness of the gravitational force, only the strongest sources can produce gravita¬ 
tional waves that are likely to be detected by present day detectors. These sources 
include colliding black holes, colliding neutron stars, and supernova explosions. The 
gravitational waves from these sources are not intrinsically weak: the metric per¬ 
turbation is of order 1 near the source. However, gravitational waves go like r~^ so 
their strength at the detector depends on the distance between detector and source. 
However black hole collisions (and neutron star collisions and supernovae) are very 
rare events; so rare that one cannot expect even one such event to occur in our galaxy 
within a reasonable amount of observing time. Thus when gravity waves are directly 
detected, the source will most likely be outside our galaxy. At such a large distance 
the r~^ falloff of the wave will insure that the amplitude of the wave at the detector 
is extremely small. The amplitude of the memory is not weaker than that of other 
aspects of the gravitational wave signal. However, there are two effects that make 
other aspects of the gravitational wave signal easier to detect: matched hltering and 
the properties of seismic noise. When black holes collide, they do so because they 
have been in orbit around each other with the orbit becoming ever smaller due to loss 
of energy from the radiation of gravitational waves, until hnally the orbit becomes 
so small that the black holes merge. (The same is true for neutron star collisions). 
Through a combination of approximation methods and numerical simulations, pre¬ 
dictions are made for the waveform for the inspiral and merging (for both the black 
hole and neutron star cases). The technique of matched hltering is to compare the 
detector signal to the predicted waveform, with a match between the two constituting 
a detection. This is a much more sensitive detection strategy than simply examining 
the detector signal, since it makes use of the large amount of information in the pre¬ 
dicted waveform. In constrast, memory is essentially only one piece of information 
and thus not subject to the effective boost in signal that matched hltering provides. 
“Seismic noise” refers to the fact that even with the best available technology the 
detector cannot be completely isolated from extraneous vibration. This vibration is 
larger at low frequencies and thus makes it more difficult to detect a low frequency 
gravitational wave. For these purposes the “frequency” of the memory signal is the 
inverse of the time that it takes the gravitational wave train to pass. Therefore mem¬ 
ory is lower frequency than the other aspects of the gravitational wave signal and 
is therefore more subject to interference by seismic noise. These issues of matched 
hltering and seismic noise are less severe for the case of supernovae than for black 
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hole or neutron star collisions. However, even in the supernova case memory is harder 
to detect than other aspects of the gravitational wave. Thus, one can expect that 
LIGO will detect gravitational waves before it detects their memory. Perhaps the 
best prospect for detecting gravitational wave memory comes from pulsar timing. EH 
Here the frequency of the gravitational waves to be measured is much lower than for 
LIGO and the sources of noise are very different. 

7 Cosmological Waves 

In the treatment of this chapter up to now we have considered memory as a prop¬ 
erty of gravitational waves in an asymptotically flat spacetime, i.e. one that ap¬ 
proaches Minkowski spacetime at large distances from the sources. However, we live 
in an expanding universe, described at the largest scale by the Friedmann-Lemaitre- 
Roberston-Walker (FLRW) metric. How much difference does that make for the 
properties of gravitational waves and their memory? Gertainly it makes a large dif- 
fence for the waves detected by BIGEP2[T]. These waves owe their existence and 
properties to the expanding universe: they start out as quantum fluctuations of the 
gravitational field that are then expanded to enormous size during a time of rapid 
exponential expansion of the very early universe. And these waves are (indirectly) 
detected through their effect on the cosmic microwave background, the light left over 
from the Big Bang explosion. But what about more ordinary gravitational waves 
and their memory? such as those produced by colliding black holes and neutron stars 
and by supernova explosions. Is the memory from these sources affected by the fact 
that they are produced in an expanding universe? One such effect is a redshift: the 
stretching of the wavelength of the gravitational wave due to the expansion of the 
universe. This effect is treated e.g. in the analysis of |32]. However, there may be 
additional effects of an expanding universe that change the properties of gravitational 
wave memory. The treatment of gravitational wave memory in an expanding universe 
is a current topic of active research of the present authors. 
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